Non-flow correlations and elliptic flow fluctuations 
in Au+Au collisions at ^/s NN = 200 GeV 
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This paper presents results on event-by-event elliptic flow fluctuations in Au+Au collisions at 
^/s NN = 200 GeV, where the contribution from non-flow correlations has been subtracted. An 
analysis method is introduced to measure non-flow correlations, relying on the assumption that 
non-flow correlations are most prominent at short ranges (|A?7| <2). Assuming that non-flow cor- 
relations are of the order that is observed in p+p collisions for long range correlations (|Ar?| >2), 
relative elliptic flow fluctuations of approximately 30-40% are observed. These results are consistent 
with predictions based on spatial fluctuations of the participating nucleons in the initial nuclear 
overlap region. It is found that the long range non-flow correlations in Au+Au collisions would have 
to be more than an order of magnitude stronger compared to the p+p data to lead to the observed 
azimuthal anisotropy fluctuations with no intrinsic elliptic flow fluctuations. 



I. INTRODUCTION 

The characterization of the collective flow of produced 
particles by their azimuthal anisotropy has proven to 
be one of the more fruitful probes of the dynamics of 
heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC). Flow is sensitive to the early stages of the 
collision and so the study of flow affords unique insights 
into the properties of the hot and dense matter that is 
produced, including information about the degree of ther- 
malization and its equation of state £Q. 

Elliptic flow, quantified by the second coefficient, «2, 
of a Fourier decomposition of the azimuthal distribution 
of observed particles relative to the event-plane angle, 
has been studied extensively in collisions at RHIC as 
a function of pseudorapidity, centrality, transverse mo- 
mentum, center-of-mass energy and system size [2H7]- A 
detailed comparison of these results to theoretical mod- 
els requires a quantitative understanding of the contri- 
butions of other many-particle correlations, referred to 
as "non-flow" and event-by-event elliptic flow fluctua- 
tions [8]. In particular, the measurement of event-by- 
event fluctuations can pose new constraints on the mod- 
els of the initial state of the collision and its subsequent 
hydrodynamic evolution [TU] . 

Comparison of the elliptic flow measurements in the 



Au+Au and Cu+Cu systems at RHIC suggests the exis- 
tence of large fluctuations in the initial geometry of heavy 
ion collisions [4]. These initial state fluctuations are ex- 
pected to lead to event-by-event fluctuations in the mea- 
sured elliptic flow signal. The measurement in Au+Au 
collisions of dynamic fluctuations in vi 1 including con- 
tributions from event-by-event elliptic flow fluctuations 
and non-flow correlations, has yielded results which are 
consistent with this expectation 

Different methods have been proposed to reduce the 
contribution of non-flow correlations to the elliptic flow 
measurements [T2l [13] . However, the application of these 
methods to the measurement of elliptic flow fluctuations 
is limited due to the complicated interplay between non- 
flow correlations and elliptic flow fluctuations [TTJl [H] . 

Ollitrault et al. have suggested estimating the magni- 
tude of non-flow from measurements of correlations in 
p+p collisions [14]. However, this estimation may not 
be completely reliable since a richer correlation structure 
is observed in Au+Au collisions at RHIC in comparison 
to the p+p system (e.g. [15] HB]). We propose a method 
to separate flow and non-flow contributions to the sec- 
ond Fourier coefficient of azimuthal particle pair distri- 
butions by studying the three-dimensional two-particle 
correlation function in (771, r]2, A</>) space. This separa- 
tion relies on the assumption that non-flow correlations 
are most prominent in short range (A?y= |r/i — r^l < 2). 
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The presumably small long range (\Ai]\ > 2) non-flow 
correlations are estimated using p+p data, and HIJING 
and PYTHIA models. Estimation of non-flow correla- 
tions using these assumptions allows the subtraction of 
the contribution of non-flow correlations to the measured 
dynamic v 2 fluctuations to obtain event-by-event elliptic 
flow fluctuations. 

This paper is organized as follows. The experimen- 
tal data is described in Section [TT| The measurement of 
the non-flow correlations and the corresponding event- 
by-event elliptic flow fluctuations are presented in Sec- 
tions HIT] and [IVJ Discussion and conclusions are included 
in Sect. |V| The numerical relation between dynamic v 2 
fluctuations, elliptic flow fluctuations and non-flow cor- 
relations is addressed in Appendix |A"} 

II. EXPERIMENTAL DATA 

The data presented here for Au+Au collisions at 
^/s NN = 200 GeV were collected during RHIC Run 4 
(2004) using the PHOBOS detector [IS]. The primary 
event trigger requires a coincidence between the Paddle 
Counters, which are two sets of sixteen scintillator de- 
tectors located at 3.2 < \rj\ < 4.5. An online vertex is 
determined from the time difference between signals in 
two sets of 10 Cerenkov counters located at 4.4 < |ry| < 
4.9, to select collisions that are close to the nominal ver- 
tex position z vtx = along the beam-axis. 

Offline vertex reconstruction makes use of information 
from different sub-detectors. Two sets of double-layered 
silicon Vertex Detectors (VTX) are located below and 
above the collision point. PHOBOS also has two Spec- 
trometer arms in the horizontal plane used for tracking 
and momentum measurement of charged particles. For 
events in the selected vertex region, the most accurate z 
(along the beam) and y (vertical, perpendicular to the 
beam) positions are obtained from the Vertex Detector, 
while the position along x (horizontal, perpendicular to 
the beam) comes primarily from the Spectrometer. 

The collision centrality is defined through bins of frac- 
tional total inelastic cross section, determined using the 
energy deposited in the Paddle Counters. In this pa- 
per, we report results for 6-45% most central events, for 
which measured dynamic v 2 fluctuations values are avail- 
able [llj . About 4 million collision events were selected 
in this centrality range by requiring that the primary col- 
lision vertex falls within \z vtx \ < 6 cm. 

The analysis presented in this paper is performed 
using the reconstructed hits in the large-acceptance 
PHOBOS Octagon silicon array, covering pseudorapid- 
ity — 3<?/ < 3 over almost the full azimuth. The angu- 
lar coordinates (77, <f>) of charged particles are measured 
using the location of the energy deposited in the single- 
layer silicon pads of the Octagon. After merging of sig- 
nals in neighboring pads, in cases where a particle trav- 
els through more than a single pad, the deposited en- 
ergy is corrected for the angle of incidence, assuming 



that the charged particle originated from the primary 
vertex. Noise and background hits are rejected by plac- 
ing a lower threshold on the corrected deposited energy. 
Depending on 77, merged hits with less than 50-60% of 
the energy loss expected for a minimum ionizing particle 
are rejected [20]. Since the multiplicity array consists of 
single-layer silicon detectors, there is no pt, charge or 
mass information available for the particles. All charged 
particles above a low-py cutoff of about 7 MeV/c at r/=3, 
and 35 MeV/c at 77=0 (which is the threshold below 
which a charged pion is stopped by the beryllium beam 
pipe) are included on equal footing. 

III. MEASUREMENT OF NON-FLOW 
CORRELATIONS 

If the only correlations between particles are due to el- 
liptic flow, then the distribution of the azimuthal angular 
separation between particles (A0 = ^i — cf> 2 ) is given by 
I + 2V r cos(2A</>), where V — v 2 {r]i) x v 2 (rj 2 ). In general, 
the second Fourier coefficient of the A<p distribution has 
contributions from both flow and non-flow correlations. 

Flow and non-flow contributions can be separated with 
a detailed study of the 77 and Ar/ dependence of the A<j) 
correlation function. Consider the distribution of A<f> be- 
tween particles selected from two 77 windows centered at 
771 and 7/2. We define the quantity uK^i!^) as the sum 
of flow and non-flow contributions to the second Fourier 
coefficient of the normalized A(f> distribution: 

V2(m,Va) = (cos(2A<£))(r>i,%) (1) 

The contributions to the second Fourier coefficient of 
the A<j) distribution can be parameterized as 

(cos(2A^)) = <7; 2 2 ) flow + ( 5, (2) 

where S is the contribution of non-flow correlations [21] . 
Using the fact that elliptic flow leads to a correlation 
between all particles in the event and creates a signal 
which only depends on pseudorapidity («2 (??)), we can 
write: 

vtivum) = v 2 (vi ) x v 2 (m) + 5(771, 772), (3) 

The measurement of non-flow correlations is therefore 
achieved in two steps, described in the following sections. 
First we measure the three dimensional (771,772, A0) cor- 
relation function to obtain v 2 (rji, 772). Then we separate 
the observed v 2 (i]i,r] 2 ) distribution to its flow and non- 
flow components. 

A. Two particle correlations analysis 

Two particle correlations have been studied extensively 
in (A77, Acf>) space using the PHOBOS detector for vari- 
ous collision systems [TCI [22- In this analysis, we extend 
the same analysis procedure to (771 , 772, Atfi) space. 
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The inclusive two-particle correlation function 
(771, 772, A(j>) space is denned as follows 



-Rn (7?1,772, A(j>) 



p mixcd (77 1 , 772,A0) 



1 



(4) 



where Pn{r]i, 772, A</>) (with unit integral in each 771,772 
bin) is the foreground pair distribution obtained by tak- 
ing two particles from the same event, then averaging 
over all pairs in all events and p mixed (77i, 772, A<fi) (with 
unit integral in each 771 , 772 bin) is the mixed-event back- 
ground distribution constructed by randomly selecting 
two particles from two different events with similar ver- 
tex position and centrality, representing a product of two 
single particle distributions. A vertex bin size of 0.2 cm 
is used in the event-mixing. 

The high occupancies measured in A+A collisions re- 
quire us to account for the high probability of multiple 
particles hitting a single pad. Furthermore, secondary ef- 
fects, such as ^-electrons, 7 conversions and weak decays, 
cannot be all rejected directly. Corrections for the high 
occupancy in the Octagon detector and the secondary ef- 
fects have been applied in the same way as in the previous 
A77, Acf> correlation analyses [TBI I22j . 

To correct for the effects of occupancy, each hit is as- 
signed a weight while calculating the correlation function. 
The weight is calculated using the centrality of the event 
and pseudorapidity of the hit (which determine the like- 
lihood of multiple particles passing through a pad for a 
given dE/dx value) and the dE/dx information. The de- 
tails of the occupancy correction can be found in Ref . [16] . 

To correct for the secondary detector effects in the 
data, correlation functions were calculated for different 
Monte Carlo event generators (PYTHIA, HIJING and 
a modified PYTHIA in which all intrinsic correlations 
have been removed) at y/s NN = 200 GeV both at the gen- 
erator level for true primary charged hadrons and with 
the full GEANT detector simulation and reconstruction 
procedure. The overall correlation structure for the re- 
constructed Monte Carlo events consists of both intrinsic 
and secondary correlations and these two sources of cor- 
relations were found to be largely independent of each 
other, i.e. the correlation from secondaries is mostly de- 
termined by sensor thickness, detector geometry, known 
cross-sections and decay kinematics [2"2"] . 

The final correlation function, R^ a ^ al (r]i, 772, A(f>) 
is calculated from the raw correlation function, 
R° a ^ w (j)x, 772, A<j>) by subtracting the contribution from 
secondary correlations: 
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data 
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(771,772, A</>) 



R 



data 



,(771,772, A4>) - 5(771,772, Acf>), (5) 



where the correction factor 5(771, 772, A(f>) is calculated by 
comparing the generator level correlation function ex- 
cluding particles outside the PHOBOS detector accep- 



tance, R 



MC 

n pri,acc 



(r)i, 772, A</>), to the correlation function 



CM CM 
> 



0.004 



0.002 




obtained with the full GEANT detector simulation and 



FIG. 1: Second Fourier coefficient of the correlation function 
R„(A(j>, 771, 772) as a function of 771 and 772 for the 40-45% cen- 
tral Au+Au collisions at ^/s NN = 200 GeV. The ridge along 
771 = 772 represents the region where non-flow contributions 
are most prominent. 



reconstruction procedure, -Rjf§ m (?7i) ?72> A0): 

5(771,772, A0) = 

R^ m (771, 772, A(f>) - l^ ri ,acc (771 , 772 , A<f>) . (6) 

The correction factor 5(771, 772, A(j>) is calculated sepa- 
rately for each centrality bin using a set of HIJING events 
with appropriate average multiplicity. More details on 
the correction factor 5 and its dependence on A77 and 
A<j) can be found in [22] . 

The second Fourier coefficient of the normalized A0 
distribution is calculated from the correlation function 
by a fit in each (771,772) bin: 

itfwfa.'Ja.Atf) - 27^(77!, 772) cos(2A</>). (7) 

The value of 7^(771,772) can also be calculated directly as 

v 2 2 (m,V2) = J flfw(M, A</>) cos(2A0)dA0. (8) 

The two methods of calculating 7^(771,772) are found to 
be equivalent within the systematic uncertainties of the 
measurement. The resulting 7*1(771,772) distribution for 
40-45% centrality bin is shown in Fig. [I] 

B. Separation of flow and non-flow contributions 

The measured 7^(771, V2) signal in Fig. [l] shows the fea- 
tures expected from Eq. [3] a ridge along A77 = where 



FIG. 2: Flow (left) and non-flow (right) components of 1)2(771, 
correlations at \Arj\ >2 are negligible. 



>T2, 



in Fig. [I] obtained by Eq. [9] and 10 assuming non-flow 



the non-flow signal is most prominent which sits on a 
plateau which can be factorized in 771 and 772. Assuming 
non-flow correlations are small at large A77 separations, 
it is possible to separate the (771, 772) *° ^ s now anc ^ 
non-flow components. 

We start by assuming that non-flow correlations at 
IA77I >2 (<5|a?7|>2) are zero. Then, we can perform a fit 



v z 2 (Vi,m) = Mvif 1 x v 2 { m ) m ■ \ m - m \ > 2, (9) 

where the fit function f2(?7) fit is an eighth order even 
polynomial. The fit in the selected A77 region can be 
used to extract the magnitude of correlations due to flow, 
^2(?7i) fit x v 2{V2) Rt , in the whole pseudorapidity accep- 
tance. Subtracting the correlations due to flow, we can 
extract the contribution of non-flow correlations: 



v 2 (m) x Mm) 



in 



(10) 



The two components of the v^rjifTfe) distribution in 
Fig. [l] are shown in Fig. [2] 

Different flow measurements with different methods 
and pseudorapidity acceptances are influenced differently 
by the non-flow correlation signal. To calculate the ef- 
fects of non-flow correlation on the measurement of dy- 
namic v 2 fluctuations performed by PHOBOS [11] , we 



calculate the average of the S(rji. 
tributions over all particle pairs: 



772) and v£{rn,r]2) dis- 



(-1) = 



(11) 
(12) 



where dN/drj is the observed charged-particle pseudora- 
pidity distribution in the PHOBOS detector. To cancel 
scale uncertainties in these quantities, we calculate the 
"non-flow ratio" given by (5) / (1*2 )■ 

The systematic uncertainty has been evaluated for the 
various stages of the non-flow ratio calculation including 
the calculation of the correlation function and the fit to 
v\ (771,772) to obtain the non-flow ratio. A "digital" occu- 
pancy correction with only the event-by-event hit density 
distribution and no dE/dx information has been used. 
Hits on the PHOBOS Vertex detector, which has a dif- 
ferent granularity from the Octagon detector have been 
added to the analysis. Monte Carlo samples with differ- 
ent average multiplicity from the data have been used in 
the correction procedure. The A77 cut used in the fit has 
been varied between 1.2 and 2.7 1 . Different fit functions 
V2(i]) iit have been used from second order up to eighth or- 
der polynomials. Finally the complete analysis chain has 
been performed by dividing the data set into 6 x 2cm wide 
vertex bins. Systematic errors are estimated for different 
steps in the analysis using the variation in the results 
with respect to the baseline due to these changes in the 
analysis. The errors in the different steps are added in 
quadrature to obtain the 90% confidence interval on the 
measurement of non-flow ratio. 

So far, we have assumed that long range (IA77I > 2) 



The Octagon detector with a pseudorapidity coverage of 
— 3<r? < 3 allows particle pairs to be studied up to Ar? = 6. 
However, in this study the Ar/ cut is constrained to At? c < 3 
such that particles from all 77 values contribute in the fit to ob- 
tain V2(rj). 
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FIG. 3: Measured value of the non-flow ratio ({5} as a function of the Arj cut (Ar/c) where non-flow correlations are 

assumed to be zero for | A.77I > Arj c for different centrality bins. The black circles (one for each panel) show values for Ar/ C = 2.1 
with the gray band denoting the 90% C.L. systematic errors on those results as described in the text. The gray squares show 
values for 1.2 < Ar] c < 2.7, which are used in the systematic error estimation. Open squares show values for Arj c < 1. 



non-flow correlations can be neglected. However, studies 
of the correlation function in p+p collisions show that 
non-flow correlations do extend out to |A?7| > 2 in ele- 
mentary collisions [22]. Furthermore, a rich correlation 
structure in high pr-triggered correlations that extend 
out to IA77I > 2 has been observed in 200 GeV Au+Au 
collision at RHIC [T5] after the estimated flow signal is 
subtracted. However, due to the inherent uncertainty 
in the flow subtraction, it is not possible to determine 
the second Fourier coefficient of this correlation struc- 
ture precisely. 

The study of the non-flow ratio as a function of the 
Arj cut (Ar] c ) for the f2(?y) fit fit carries important infor- 
mation on the magnitude of non-flow at large Aij sepa- 
rations. If non-flow correlations are short ranged, we ex- 
pect that the fits should yield non-flow ratio results that 
saturate for large values of Ar] c . The extracted value of 
(5) I (vf ) is plotted as a function of the A?y c , where it is 
assumed that 5 is zero for |A7y| > Arj c , for different cen- 
trality bins in Fig. [3] The saturation expected if non-flow 
correlations are short-range is indeed observed. However, 
it should be noted that the same saturation pattern could 
also be observed with a finite magnitude of non-flow that 
has little Aij dependence in the region A77 > 1.2. 

To quantitatively assess the effect of non-zero non-flow 
correlations at large Arj separations, we analyze the cor- 
relation functions obtained from Monte Carlo event gen- 
erators. In p+p collisions, the magnitude of non-flow 
correlations, S, can be directly calculated as the second 
Fourier coefficient of A0 correlations since elliptic flow is 
not present [22]. If A+A collisions were a superposition 
of p+p collisions, the value of S would be diluted due to 
the presence of uncorrelated particles. To compare the 
strength of non-flow correlations in HIJING (Au+Au) 



and PYTHIA (p+p) models and p+p collisions, we cal- 
culate the value of 5 scaled by the average event multi- 
plicity, shown in Fig. [4p. Both models are observed to 
roughly reproduce the strength of non-flow correlations 
in p+p collisions at large Ar/. Due to large systematic 
uncertainties in the p+p data, HIJING simulations are 
used to model the long range non-flow correlations in 
Au+Au collisions by assuming non-flow correlations in 
data are some multiplicative factor, to, times the non- 
flow in HIJING (S M c(Vi,V2)) for |A?y| > 2. This can be 
incorporated by modifying Eq. [9j 

"I {Vl > ^2) - mSuc (Vi , m ) = 

v 2 (r ]1 f t x v 2 { m f t ■ \Ar,\>2. (13) 

The resulting non-flow ratio, (5) / (wf), found by ap- 
plying Eqs. |l0l|T2"1 with the modified 1)2(77) results, is 
plotted as a function of centrality in Fig. [5] for different 
assumptions on the magnitude of non-flow at | Ai]\ > 2. If 
non-flow correlations are assumed to be present only in 
IA77I < 2 (m = 0), it is found that they account for ap- 
proximately 10% of the observed v\ signal averaged over 
1 77 1 < 3. The results do not change significantly if the 
long range non-flow correlations (<5|Ai;|>2) are taken to be 
the same as the correlations in HIJING (to = 1 instead 
of to = 0). 

The upper limit on the non-flow ratio, also shown in 
Fig. [5] is drawn from the measurement of dynamic v 2 



2 The large uncertainty in the p+p data at Ar; = is due to 
<5-electrons and 7 conversions, which may not be completely de- 
scribed by GEANT simulations [22]. 
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FIG. 4: The magnitude of non-flow correlations (8) scaled 
by the charged particle multiplicity (n) in the pseudorapid- 
ity range \rj\ < 3 as a function of particle pair pseudorapid- 
ity separations (Arj) for p+p data and different Monte Carlo 
generators with no flow correlations at ^/s NN = 200 GeV. The 
results for p+p data (squares) with 90% C.L. systematic er- 
rors are obtained from two particle Ar], A(f> correlations [22] . 
Statistical errors are not shown. 

fluctuations [TT] assuming that the observed fluctuations 
are all due to non-flow correlations. The calculation of 
this limit is described in Appendix [X] This limit corre- 
sponds to non-flow correlations in Au+Au collisions that 
are more than an order of magnitude higher than the 
expected correlations from p+p collisions for \Ar)\ > 2 
(to > 10). 

IV. ELLIPTIC FLOW FLUCTUATIONS 

An event-by-event measurement of the anisotropy in 
heavy ion collisions yields fluctuations from three sources: 
statistical fluctuations due to the finite number of parti- 
cles observed, elliptic flow fluctuations and non-flow cor- 
relations. We have previously measured the dynamic 
fluctuations in v 2 by taking out the statistical fluctua- 
tions with a study of the measurement response to the 
input v 2 signal [11] . The new results on the magnitude 
of non-flow correlations presented in the previous section 
can be used to decouple the contributions of genuine el- 
liptic flow fluctuations and non-flow correlations to the 
measured dynamic fluctuations. 

Let us denote the observed distribution of the event- 
by-event anisotropy as g(v 2 hs ), the distribution of the 
intrinsic elliptic flow value as f{v 2 ) and the expected dis- 
tribution of v 2 hs f° r a Axed value of v 2 as K(v 2 hs , v 2 ). We 
assume f{v 2 ) to be a Gaussian in the range v 2 > with 
two parameters, mean ((v 2 )) and standard deviation (tr). 
The dynamic fluctuations in v 2 , can be calculated by un- 
folding the experimental measurement g exp (v2 hs ) with a 
response function K^ p (v 2 hs , v 2 ) which accounts for de- 



F1G. 5: The non-flow ratio «5) / («f)) in the PHO- 
BOS Octagon detector acceptance as a function of num- 
ber of participating nucleons (Af part ) in Au+Au collisions at 
^/Snn = 200 GeV. The black squares show the results with 
the assumption that non-flow correlations are negligible at 
IA77I > 2. The shaded band shows the 90% confidence sys- 
tematic errors. The lines show different assumptions about 
non-flow at [A77I > 2. The open circles with 90% C.L. sys- 
tematic errors, show the upper limit on {8} / (vf) obtained by 
assuming that the measured dynamic fluctuations in V2 are 
due to non-flow alone. 

tector effects and statistical fluctuations: 

5 exp (u ob S) = I ^exp (u ob S) y 2 ) f dyn{v2 ) dv2 . {U) 
JO 

The calculation of intrinsic flow fluctuations (/fl ow (^2)) 
from measured dynamic fluctuations {fdyn(v 2 )) can be 
summarized by the following equation: 

/ K n (v° hS ,V2)f d yn(v 2 )dv 2 

Jo 

= [ K n , s (v% hs ,v 2 )f Row (v 2 )dv 2 , (15) 
Jo 

where K n (v 2 hs ,v 2 ) and K ny s(v 2 hs , v 2 ) are the response 
functions for an ideal detector with and without non-flow 
correlations respectively. Equation [15] gives the distribu- 
tion of observed anisotropy for an ideal detector g{v 2 hs ), 
such that on the left hand side the non-flow correla- 
tions are encoded in the dynamic v 2 fluctuations, and 
on the right hand side, they are accounted for in the re- 
sponse function K n . s {vf s , v 2 ). The response functions 
K n (v 2 bs ,v 2 ) and K ny s(v 2 hs , v 2 ) are given by a Bessel- 
Gaussian distribution |23| defined as 

BG {vf s -v 2 ,a s ) = ^ 

xex H - Li Jn^rJ' (16) 
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flow fluctuations are found to be roughly 30-40% if the 
magnitude of non-flow correlations are assumed to be 
small for | A77I > 2. The observed values of relative elliptic 
flow fluctuations correspond to 87-97% (79-95%) of the 
previously measured dynamic v% fluctuations if non- 
flow correlations at \Arj\ > 2 are assumed to be zero (three 
times the magnitude in HIJING). 

Also shown in Fig.[6]are relative fluctuations in the par- 
ticipant eccentricity obtained from MC Glauber [11] and 
color glass condensate(CGC) [21] calculations. The mea- 
sured values of elliptic flow fluctuations are observed to 
be consistent with both models over the centrality range 
under study if the long range non-flow correlations are 
neglected. The same conclusion holds if the long range 
correlations are assumed to be three times stronger than 
in p+p collisions, as modeled by HIJING. 



FIG. 6: Relative elliptic flow fluctuations (ofl ow / (v2) aow ) as 
a function of number of participating nucleons (iV par t) in 
Au+Au collisions at ^/s NN = 200 GeV. The black circles show 
the results with the assumption that non-flow correlations are 
negligible at \Arj\ > 2. The shaded band shows the 90% con- 
fidence systematic errors. The thin lines show results for dif- 
ferent assumptions on the magnitude of non-flow at \Arj\ >2. 
The continuous and dashed thick lines show a(e par t) /{e pa rt) 
values calculated in Glauber MC [TT] and CGC [24] models, 
respectively. 



where Iq is the modified Bessel function. The fluctuation 
term a s in the response function is a quadratic sum of sta- 
tistical fluctuations (a n — \j\p2m) due to finite number 
of particles (n) observed in the detector and a contribu- 
tion from non-flow correlations (as = y/S/2). 

Equation 15 cannot be simplified analytically. How- 
ever, it can be solved numerically to calculate relative 
elliptic flow fluctuations (er now / («2)fl ow ) that correspond 



to the measured dynamic vi fluctuations (odyn/ (^2)) and 
the non-flow ratio ((5) / (v%)) for different assumptions 
on non-flow at |Ary| > 2. The details of the numerical 
calculation are given in Appendix [S] It has been sug- 
gested that the relation between these quantities can be 
approximated as (Tj yn = + ct| ow [TJ] . We have found 
that this approximation does not hold in the range of our 
experimental results (<Tdyn/ (^2) > 0.3). 

The systematic error in the magnitude of relative ellip- 
tic flow fluctuations is obtained by propagating the errors 
in the measured quantities fXdyn 

/ (v 2 ) and (S) / (vf) an d 
by varying the procedure to calculate cx now / («2)fl ow from 
these quantities. The errors from different sources are 
added in quadrature to obtain the 90% confidence in- 
terval. The error propagated from the uncertainty in 
Cdyn/ ( v 2) is the dominant contribution to the uncer- 
tainty in CT flow / (v 2 ) flow - 

The relative fluctuations in the event-by-event elliptic 
flow, corrected for contribution of non-flow correlations 
are presented in Fig. [6] as a function of the number of 
participating nucleons, in Au+Au collisions at ^/s NN = 
200 GeV for 6-45% most central events. The elliptic 



V. SUMMARY AND CONCLUSIONS 

We have presented new data on the magnitude of non- 
flow correlations and the event-by-event elliptic flow fluc- 
tuations corrected for non-flow correlations in Au+Au 
collisions at ^/s^ = 200 GeV. The measurement of non- 
flow correlations is achieved by utilizing a new correlation 
analysis with the assumption that non-flow correlations 
are of the order that is observed in p+p collisions for 
long range correlations > 2). The non-flow corre- 

lations averaged over the PHOBOS Octagon acceptance 
(— 3 < 77 < 3) are found to be large, constituting approxi- 
mately 10% of the measured v\ signal. Studying the de- 
pendence of expected azimuthal anisotropy fluctuations 
due to non-flow correlations, it is found that the long 
range non-flow correlations in Au+Au collisions would 
have to be more than an order of magnitude stronger 
compared to the p+p data for non-flow correlations to 
lead to the observed azimuthal anisotropy fluctuations 
with no intrinsic elliptic flow fluctuations. The method 
presented in this paper can be generally applied in large 
acceptance detectors to study the contribution of non- 
flow correlations to the flow signal measured with differ- 
ent approaches. 

The magnitude of event-by-event elliptic flow fluctua- 
tions were calculated by subtracting the contribution of 
non-flow correlations to the measured values of dynamic 
V2 fluctuations. If the inclusive long range non-flow cor- 
relations in A+A collisions are assumed to be of the or- 
der of magnitude that is observed in p+p collisions, the 
magnitude of event-by-event elliptic flow fluctuations are 
found to be in agreement with predicted fluctuations of 
the initial shape of the collision region in both Glauber 
and Color Glass Condensate models. Therefore these re- 
sults support conclusions from previous studies on the 
importance of geometric fluctuations of the initial col- 
lision region postulated to relate elliptic flow measure- 
ments in the Cu+Cu and Au+Au systems [1]. 

This work was partially supported by U.S. DOE 
grants DE-AC02-98CH10886, DE-FG02-93ER40802, 
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Appendix A: Numerical calculations relating 
measured quantities to elliptic flow fluctuations 

In this section, we describe the numerical calculations 
performed to relate the measured values of dynamic v 2 
fluctuations (cdyn/ (^2)) and non-flow ratio ((S) / (v%)) 
to intrinsic elliptic flow fluctuations (crfl ow / ( t '2)fl ow )- 

We start by assuming the mean value of the ellip- 
tic flow distribution and the magnitude of statistical 
fluctuations to be given as («2)fi ow = 0.06 and a n = 
0.6 x («2)fl ow = 0.036 (see Eq. 16 1. Then, for given values 
// ( w 2)fl ow and (5) / (v 2 ), the expected distribution 



of cr flo 

of the observed event-by-event anisotropy t>2 bs can be 
calculated as 



<?K bs ) 



(Al) 



where fa ow {v 2 ) is a Gaussian in the range v 2 > with 
mean and standard deviation values given by (v 2 ) Row and 
% w , respectively, and K n j(v% hs , v 2 ) is given by a Bessel- 
Gaussian (see Eq. 16), 



K n , s (v° hs ,v 2 )=BG (v° bs ;v 2 ,a s ) 



(A2) 



The fluctuations encoded in the response function 
K n ,s{v 2 hs , v 2 ) are given as a 2 s = a\ + <r|, where as can 
be calculated from (w2) flow , <7fk>w and (S) / (v 2 ): 



2aj = (S) 



(S)x 
(S)/(v 2 2 ) 

l-W/(«2 2 > 



(A3) 
(A4) 



x(WL + ^ow)- (A5) 



In this derivation, it has been noted that the (v 2 ) defined 
in Eq. 12 includes contributions from flow fluctuations 



and non-flow correlations. 

Next, we calculate the dynamic fluctuations in the 
measured ^2 bs distribution, g{v'^ os ) 1 by using a response 
function which incorporates only statistical fluctuations 
but not non-flow correlations, 



sK bs ) = 



K n (v% hs , v 2 )f dyn (v 2 )dv 2 . 



(A6) 



Assuming the dynamic v 2 fluctuations are described by 
a Gaussian, fdyn{v 2 ), in the range v 2 > with mean and 
standard deviation values given by (v 2 ) and cr^ ym t he 
value of Cdyn/ (^2) can be obtained by fitting Eq. |A6| 

The resulting distribution of 0d yn / (^2) as a function 
of cr flow / (w 2 ) flow and (S) / (vf) is shown in Fig. [7] The 
value of crflow/ ( v 2)a ow corresponding to measuredvalues 
of (Tdyn/ {v 2 ) and (S) / can be extracted from this 
distribution. Furthermore, the values for <7fl ow / ( v 2)f{ ow — 
can be used to set an upper limit on the magnitude of 
the non-flow ratio. 

Since the related quantities are given as ratios, the 
value of (t>2)fl ow set at the beginning is arbitrary. It was 
observed that a n / (v 2 ) is roughly given by 0.6 for the 
dynamic v 2 fluctuations measurement for all centrality 
bins in the centrality range studied. The calculation was 
repeated for values of a n j (f2)fl ow =0.4 and 0.8. The dif- 
ferences in results, which were found to be small, are 
incorporated in the systematic errors. 



[1] P. F. Kolb, P. Huovinen, U. W. Heinz, and H. Heiselberg, 

Phys. Lett. B500, 232 (2001). 
[2] B. B. Backet al. (PHOBOS), Phys. Rev. Lett. 94, 122303 

(2005). 

[3] B. B. Back et al. (PHOBOS), Phys. Rev. C72, 051901 
(2005). 

[4] B. Alver et al. (PHOBOS), Phys. Rev. Lett. 98, 242302 



(2007). 

[5] B. B. Back et al. (PHOBOS), Nucl. Phys. A757, 28 
(2005). 

[6] J. Adams et al. (STAR), Nucl. Phys. A757, 102 (2005). 
[7] K. Adcox et al. (PHENIX), Nucl. Phys. A757, 184 
(2005). 

[8] H. Song and U. W. Heinz, J. Phys. G36, 064033 (2009). 



9 



[9] T. Osada, C. E. Aguiar, Y. Hama, and T. Kodama 
(2001), arXiv:nucl-th/0102011vl. 
[10] B. Alver et al. (PHOBOS), Phys. Rev. C77, 014906 
(2008). 

[11] B. Alver et al. (PHOBOS), to appear in Phys. Rev. Lett 

(2010), arXiv:nucl-ex/0702036v2. 
[12] N. Borghini, P. M. Dinh, and J.-Y. Ollitrault, Phys. Rev. 

C64, 054901 (2001). 
[13] A. Bilandzic, N. van der Kolk, J.-Y. Ollitrault, and 

R. Snellings (2008), arXiv:0801.3915vl. 
[14] J.-Y. Ollitrault, A. M. Poskanzer, and S. A. Voloshin, 

Phys. Rev. C80, 014904 (2009). 
[15] J. Adams et al. (STAR), Phys. Rev. C73, 064907 (2006). 
[16] B. Alver et al. (PHOBOS), Phys. Rev. C81, 024904 

(2010). 

[17] B. I. Abelev et al. (STAR), Phys. Rev. C80, 064912 



(2009) . 

[18] B. Alver et al. (PHOBOS), Phys. Rev. Lett. 104, 062301 

(2010) . 

[19] B. B. Back et al. (PHOBOS), Nucl. Instrum. Meth. 

A499, 603 (2003). 
[20] B. B. Backet al. (PHOBOS), Phys. Rev. Lett. 87, 102303 

(2001). 

[21] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C58, 
1671 (1998). 

[22] B. Alver et al. (PHOBOS), Phys. Rev. C75, 054913 
(2007). 

[23] J.-Y. Ollitrault, Phys. Rev. D46, 229 (1992). 
[24] H.-J. Drescher and Y. Nara, Phys. Rev. C76, 041903 
(2007). 



